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ABSTRACT 



Here we report on the detailed analysis of the 7-ray light curve of a luminous 
blazar PKS 1510—089 observed in the GeV range with the Large Area Telescope 
(LAT) onboard the Fermi satellite during the period 2011 September - Decem- 
ber. By investigating the properties of the detected three major flares with the 
shortest possible time binning allowed by the photon statistics, we find a variety 
of temporal characteristics and variability patterns. This includes a clearly asym- 
metric profile (with a faster flux rise and a slower decay) of the flare resolved on 
sub-daily timescales, a superposition of many short uncorrelated flaring events 
forming the apparently coherent longer-duration outburst, and a huge single iso- 
lated outburst unresolved down to the timescale of three-hours. In the latter 
case we estimate the corresponding 7-ray flux doubling timescale to be below 
one hour, which is extreme and never previously reported for any active galaxy 
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in the GeV range. The other unique finding is that the total power released 
during the studied rapid and high-amplitude fiares constitute the bulk of the 
power radiatively dissipated in the source, and a significant fraction of the total 
kinetic luminosity of the underlying relativistic outfiow. Our analysis allows us 
to access directly the characteristic timescales involved in shaping the energy 
dissipation processes in the source, and to provide constraints on the location 
and the structure of the blazar emission zone in PKS 1510—089. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal 
— galaxies: active — galaxies: jets — quasars: individual (PKS 1510—089) — 
gamma rays: galaxies 



1. Introduction 



Blazars constitute a class of radio-loud active galactic nuclei for which the observed 
broad-band spectra are dominated b y the beamed and high ly variable emission of the in- 
nermost parts of relativistic jets (see lUrry fc Padovanilll995l ). Detailed and systematic in- 
vestigation of the blazar variability at different wavelengths provides crucial information on 
the location and the structure of the energy dissipation zone(s) in blazar jets, as well as 
on the radiative and particle acceleration mechanisms involved in the production of the ob- 
served radiation. Recently, thanks to the operation of the modern high-energy instruments 
including ground-based Cherenkov telescopes and space-borne satellites, such studies can be 
performed also in the 7-ray regime. The observational results obtained in this way so far are 
rather striking. 

Nearby low-luminosity blazars of the BL Lac type were obse rved to be variable in 
the TeV range at the extremely short timescales of a few minutes (lAharonian et al.l 120071 : 
Albert et al.l 120071 ). challenging the homogeneous synchrotron-self-Compton emission mod- 



els typically discussed in the context of the v ery high-energy (VHE ; photon energies > 



0.1 TeV) 7-ray emission of BL Lacs (see, e.g., iBegelman et al.ll2008l : iGiannios et al.ll2009l ). 
In high-luminosity blazars of the 'Flat Spectrum Radio Quasar' type (FSRQs), equally spec- 
tacular high-amplitude flux changes on the sub-daily timescales in the GeV range have been 
established during the flaring states of several sources based on the data provided by the 
Fermi Gamma- Ray Space Telescope (see §4 below). These flndings, es pecially when com- 



bined with the une xpected detections of a few FSRQs in the TeV range ( lAlbert et al.l 12008 



Aleksic et al.l 1201 if ), again put into question the validity of the one-zone inverse-Compton 
(IC) scenarios widely applied in the modeling of the high-energy (HE; g-^ — 0.1 — 100 GeV) 



7-ray continua of such objects (e.g., iBottcher et al.ll2009t iTavecchio et al.l 1201 ll ). 
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One of the FSRQs for which rapid HE outbursts have been seen with Fermi is PKS 
1510—089. This is a luminous blazar located at the redshift of z = 0.361 (luminosity distance 
c^L — 1-91 Gpc)0. The inner structure of the source is characterized by the drar natic morpho- 



logical changes with the apparent superluminal velocities up to /3app — 20 (e.g.. lOrienti et al. 



201ll ). and the extreme mis-al ignment between the milli-arcsec and arcsec-scale jets (~ 180° 
projected; iHoman et al.ll2002l ). implying ultrarelativistic velocities of the emitting plasma 
(jet bulk Lorentz factor Tj ~ 20) at small inclinations to the line of sight (jet viewing angle 
6j ^ 3°). The synchrotr on emission of the inn er jet in PKS 1510— 089 peaks around the 



infrared frequencies (e.g., iNalewajko et al.ll2012l ): the UV segment of the source spectrum is 



contr ibuted by accretion disk (most clearly during the quiescent states; iD'Ammando et al. 



20091 ). and the X-ray continuum is dominated by the IC emission of the low-energy jet 



electrons (e.g.. iKataoka et al.ll2008l ). PKS 1510— 089 is also an established 7-ray emitter, 



detected i n the HE range by the EGRET instrument on board CGRO, Ferm i, and A GILE 
satellites ( Hartman et al. 1999; Abdo et al. 2010c; D'Ammando et aD 2011 ) , as well as i n 
the VHE range by the H.E.S.S. and MAGIC telescopes dWagner et aJboiol : ICortinall2012h . 

Here we report on the detailed analysis of the spectacular 7-ray outbursts of PKS 
1510—089 detected in the GeV range with LAT onboard the Fermi satellite during the period 
2011 September - December (§ 2). We study the profiles of the HE flares of PKS 1510-089 
with different time binning, finding a variety of the temporal characteristics and variability 
patterns with the flux doubling timescales as short as below one hour (§3). We discuss 
the energetics of the flares, concluding on the implication of the obtained results for the 
general understanding of the location and the structure of the energy dissipation zone in the 
relativistic jet of PKS 1510—089 and other luminous blazars of the FSRQ type (§4). 



2. Data Analysis 



The Ferm^-LAT is a pair-conversion 7-ray telescope sensitive to photon energies from 
20MeV up to 300 GeV, characterized by the energy resolution typically ~ 10% and the 
angular resolution (68% containment angle) better than 1° at energies above 1 GeV. With 
its large field of view of 2.4 sr it ope rates mostly in a su rvey mode scanning the entire sky 
every 3 hours (for further details see lAtwood et al.ll2009l ). 



The data discussed in this paper were collected from MJD 55834 (2011 September 30) 
to MJD 55903 (2011 December 8), the period which overlapped with three major 7-ray 
outbursts of PKS 1510—089. Only the events with energies greater than 100 MeV and zenith 



^We assume ACDM cosmology with Oa = 0.73, = 0.27, and Ha = 71kms~i Mpc"^ 
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angles < 100° were included in this analysis. We selected good time intervals (GTIs) by 
using a logical filter of "DATA_QUAL==1", "LAT_C0NFIG==1" , and "ABS(R0CK_ANGLE)<52". 
The LAT data collected in that manner were analyzed with an unbinned likelihood analysis 
method using the standard analysis tool gtlike, which is a part of the Fermi-LAT Science 
Tools software package (v9r27pl). The P7S0URCE_V6 set of instrument response functions 
was utilized. 

The 7-ray photons were extracted from a circular region of interest (ROI) with ra- 
dius 10°, centered at the radio position of the source. The Galactic diffuse emission tem- 
plate version "gal_2yearp7v6_v0 . f its" and the isotropic diffuse emission template ver- 
sion "iso_p7v6source.txt" were assumed in the modeling. The source model consisted 
of PKS 1510—089 and other point sources within the RO I and the surroun ding 5°-wide 
annulus taken from the second-year LAT catalog (2FGL; iNolan et al.l 120121 ). Moreover, 
we included an additional variable 7-ray point source located at (R.A., Dec.) = (233.168, 

— 13.311), i.e. about 6.4° awa y from PKS 1510—08 9. This object, tentatively associated 

with the FSRQ TXS 1530-131 (IGasparrini &: Cutinill201ll ). did not appear in the 2FGL, al- 
though it was clearly detected during the period considered here at relatively high flux level 
of F>iooMeV ^ (1.1 ±0.1) X 10-7phcm-2s-i. 

In the following analysis of the source light curves, we flxed the fluxes of the diffuse 
emission components at the values obtained by fltting the data collected over the entire 
period discussed here. Normalizations of all the point sources within 10° radius around 
PKS 1510—089 were set free, while the other relevant parameters of the neighboring 7-ray 
emitters were flxed following the second-year LAT cata log. For each time bin analyzed, 
point sources with the test statistic (TS) values < (see lMattox et al.lll996l ) were removed 
from the source model. We approximated the 7-ray continuum of PKS 1510—089 with a 
simple power-law model, keeping photon indices free during the likelihood analysis of the 
source spectra for all the time binnings considered (down to 3h). The power-law flts were 
acceptable in each case, indicating that more complex spectral sh apes, with more degrees of 
freedom, were formally not required (but see lOrienti et al.ll2012l ). 



3. Results 



Figured] presents the daily 7-ray light curve of PKS 1510— 089 at photon energies 0.1- 
300 GeV during the period analyzed in this paper. In our analysis, 95% confldence level 
flux upper limits correspond to the detection signiflcance values TS < 10, which is a con- 
ventional choice in th e analysis of daily-bi nned light curves of bright LAT sources (see, e.g.. 



Tavecchio et al.l l20ld : lOrienti et al.ll2012[ ). As shown, during the discussed time interval 
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three major high-amplitude 7-ray outbursts of the source were detected with photon fluxes 
-^>iooMcV ^ 10~^ phcm~^ s~^ (see Tabled]), and flux doubling timescales less than a day. The 
excellent photon statistics allowed us to study these flares with shorter time binning, down 
to the minimum 3h dictated by the survey mode of the LAT instrument. 

Figure|2] presents the light curves of PKS 1510—089 around the time of the first major 
7-ray outburst, binned in the intervals of 12 h, 6 h, and 3 h (upper, middle, and lower panels, 
respectively). As shown, the rising segment of the flare is unresolved down to the timescale of 
3 h, and this is a truly unique finding for a luminous blazar. Previously all the high-amplitude 
flux changes of FSRQs detected in the GeV range were characterized by longer (> Id) 
timescales, while any shorter variability consisted of a small-amplitude flickering only (see 
§4 below). Here, instead, the recorded flux increases from about Fi ^ 7 x 10~^phcm~^s~^ 
up to F2 ~ 45 X 10~^ ph cm~^ s~^ within At = 3 h, giving formally the flux doubling timescale 
of Td = At X In 2/ ln(F2/-Fi) :^ 1 h only, or equivalently the exponential growth timescale 
Td/ ln2 :^ 1.5 h. This value should be considered as an upper limit only, because of a limited 
exposure of PKS 1510—089 during the analyzed 3 h visibility window. Interestingly, the decay 
segment of the flare seems to be marginally resolved with the 3 h binning, implying the e- 
folding decay timescale of about 4h. The evaluated photon index 2.0, together wit h 



the large Compton dominance established for PKS 1510— 089 (e.g., iNalewajko et al.ll2012l ). 
imply that the bulk of the radiative energy released during the flare is contained within the 
HE range. 

The second 7-ray outburst of PKS 1510—089 for which the LAT light curves in 12, 6, 
and 3h bins are presented in Figure[3l constitutes a very different case. Here the flare seems 
to be resolved with 12 h-binning, displaying shorter exponential growth and a slower linear 
decay, as expected in most of the models of FSRQs' variability involving a fast injection of 
accelerated electrons and their slower radiative coolin g dominated by th e Comptonization of 



the soft photons produced externally to the jet (e.g., ISikora et al.ll200ll ). However, with the 
minimum 3h binning a signiflcant sub-structure of a flare becomes prominent, consisting of 
several apparently chaotic and unresolved yet still large-amplitude events, often characterized 
by the flux doubling timescales < 3 h. This clearly illustrates the fact that with the limited 
time resolution, the apparent proflles of high-energy outbursts in blazar sources may not 
reflect the exact temporal characteristics of the source flux changes. 

Finally, FigureH] presents the light curves of PKS 1510—089 around the time of the third 
major 7-ray outburst, binned again in the intervals of 12 h, 6h, and 3h. As shown, the flare 
seems to be nicely resolved in short binning, displaying a moderately asymmetric proflle 
with a faster flux increase (doubling timescale between 3h and 6h), and a longer flux decay 
(e-folding timescale of about llh). However, we cannot exclude a possibility that with even 
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shorter binning of the hght curve, this smooth and seemingly coherent flaring event would 
be decomposed into a series of rapid overlapping but not necessarily related sub-events. 



Discussion and Conclusions 



Strong HE fl ares from PKS 1510-089 have b een detected with AGILE in 2008 March 
and 2009 March fjP'Ammando et al.ll2009l . l201ll ). with the daily-integrated peak fluxes of 
F>iooMcV — (3 — 7) X 10^^ ph cm~^ s~\ photon indices ~ 2.0, and the flux doubling 
timescales of the order of a day. Similar results were reported by the F ermz-LAT Collab- 
oration for the period 2008 September - 2009 June (lAbdo et al.ll2010d ). consisting of the 
detection of several outbursts lasting for a few/several days with the daily-integrated peak 
fluxes -F>iooMcV — (2 — 8) X 10~^phcm~^ s~^, photon indices < 2.5, and the characteris- 
tic e-folding flux variability timescale estimated to be of the order of 3 h (for the minimum 
6 h-binning of the light curve applied) . 

The 2009 March and April Fermz-LAT data for PKS 1510—089 were analyzed also by 
Tavecchio et al.l (120101 ). who found significant flux changes by a factor of two or more oc- 
curring on the timescale of 6 h, with approximately symmetrical flare profiles, and also 
well-defined events characterized by even shorter variability, which however could not be 
claimed at a high confidence level due to large errors related with lower flux level of the 
source. Similar behavior with the flux doubling timescales of the order of several hours 



has been reported for some other flar ing FSRQs observe d with LA 



jAbdo et all 120091). PKS 15 02+1 06 JAbdo et all l201oi) . 3C273 



3C 454.3 (Ackermann et al. 2010: Tavecchio et al. 



20101 ). Finally, 



, namely PKS 1454 -354 
Obi), and 



Abdo et al. 



20 



Foschini et al.l ( 1201 ll ) set 



the upper limits of 2 h on the observed 7-ray doubling time scale in quasars 3C 454.3, 3C 273, 
and PKS 1222+216 during their bright flaring states with daily photon fluxes -F>iooMeV ex- 
ceeding 10~^phcm" 



-2s-^ 



The isotropic daily-averaged HE luminosity of the first flare analyzed here is 
7 X lO^^ergs"^. The corresponding total power emitted in 7-rays (i.e., the po wer as would 



be m easured by the detector completely surrounding the emitting region; e.g., ISikora et al. 



19971 ) is therefore L 



7, em 



-^^7,iso/4r? 5 X lO^'^ergs"^ (assuming Fj ~ 20; see §1), which 
is almost exactly the same as the total kinetic power of the PKS 1510—089 jet emerging 
from broad-band modeling based on different datasets and model assumptions, Lj > 5 x 
]^q45 e];p-g-i^ a]2d also as the observed UV disk luminosity in the sys tem, Ldisc — 5xl0'^^ergs~^ 
(IKataoka et al.lbood : IP'Ammando et al.lbood : kbdo et al.lboiOd ). This implies that, during 
the discussed flaring event, the power dissipated in the jet within less than a day and emitted 
as 7-ray photons constitutes the bulk of the total kinetic luminosity carried out by the 
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j^,cm/L^ < 1, and also a substantial fraction of the entire available accretion power, 
em/ Lace — 0.1 (assuming the standard ?7disk — 10% radiative efficiency for the accretion 
disk, L^r 



outflow, L 

Ldisc/'7disk — 5 X 10^6 e rg s~^). Note in thi s context that, for the black hole mass 
in the system A^bh =^ 5 x lO^M© jAbdo et al.lboiOch . the active nucleus in PKS 1510-089 
accretes at the maximum Eddington rate. Lace ~ ^Edd- A very similar set of the source 
parameters, implying the extremely efficient c onversion of the acc retion power to the jet 7- 
ray luminosity has been established before by iTanaka et al.l (120111 ) for the analogous blazar 
PKS 1222+216 observed with LAT during its flaring state. 

Equally interesting is the analysis of the emerging timescales and the related (via the 
causality arguments) emission zone spatial scales. In particular, the observed flux doubling 
timescale of — 1 h and the bulk Lorentz factor Fj 2^ 20 (equal by assumption to the jet 
Doppler factor, consistently with the expected jet inclination 6j ~ 3°) give the spatial scale 
of the emitting region R^ar < crdrj/(l + 2) ^ 1.5 x 10^^ cm. Meanwhile, the gravitational 
radius of the PKS 1510— 089 supermassive black hole is = GA^bh/c^ ^ 7 x 10^'^ cm. 
Assuming a very basic scenario in which the scale of the event horizon sets a lower limit on 
the spatial scale of the jet disturbances that can be identifled with the zones of the enhanced 
energy dissipation, one should expect such structures, when created near the black hole, to 
be advected along the outflow and to release the bulk of their power around rem — Tj^'^g — 
3 X 10^^ cm distances from the core (IBegelman et al.ll2008l ). The characteristic radial scale of 
the outflow at that point is then expected to be approximately Rj ~ rcm/rj — 1.5 x 10^^ cm, 
following the standard expectation for the jet opening angle ~ l/Pj. The agreement between 
the derived values of -Rvar and Rj is striking. 

This identified blazar zone would be located inside the region of the highest ioniza- 
tion of the broad line-emitting circumnuclear clouds ('broad line region '; BLR), for which 
the characteristic scal e in the discussed system is tblr — 2 x 10^^ cm (lAbdo et al.ll2010c : 
Nalewajko et al.ll2012l ). There the energy density provided by the line-emitting clouds should 
exceed energy densities of the other photon fields in the jet rest frame, and hence the dom- 
inant production of the 7-ray photons should be related to the IC upscattering of the UV 
emission (obse rved energies £n — 10 eV) reprocessed with the £rt,r = 10% efficiency within 



the BLR (e.g.. iGhisellini fc Tavecchid l2009l : ISikora et al.ll2009l ). The corresponding cooling 
timescale for the electrons emitting 7-rays with the energies of = 100 MeV, as measured 
in the observer frame, would then be r^ad — (3?Tiec/4(TTMBLR) ^ [£oi'^ + /^-yY^"^ ^ lOmin, for 



the jet comoving BLR photon energy density Mblr — ^blr Ldiskrf /4vrr 



BLR*- 



10 erg cm 



This timescale is shorter by a factor of 10 — 50 than the observed e-folding decay timescales 
of the flares, implying that the observed flux decrease is shaped not solely by the radiative 
energy losses, but instead by a combination of different factors. These other factors may be 



related either to the geometry and sub-structure of the emitting region (e.g.. iTanihata et al. 
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200ll ). or to a residual particle acceleration still ongoing after the peak of a flare. 



The 'near-dissipation zone' scenario, with the dominant emission region located rela- 
tively close to the central engine (< 0.1 pc), was advocated in the literature for FS RQs in 
general based on the modeling of their HE 7-ray spectra (jPoutanen &: SternI 120101 ). The 
complication arises, however, due to the aforementioned detection of a few FSRQs, includ- 



ing PKS 1510-089 an d PKS 1222+216, at TeV photon energies (see iTanaka et al.l I2OII 



Tavecchio et al.ll201ll ). The emerging agreement is that such VHE emission, if detected, 
must be produced instead at further distances from the core, i.e. beyond the characteristic 
scale of the circumnuclear dust (> 0.1 pc). 

The broad-band variability studies enabled by the multiwavelength campaigns carried 
out during the recent years revealed rather complex behavior of PKS 1510—089 with no 
HE/X-ray correlations, weak HE/UV correlations, and pronounced HE/optical correlations; 
some of the HE/optical correlations were found to consist of optical flares lagging the 
HE outbursts by several days, while other events showed a time -lag consistent with zero 
(IMarscher et al.ll2010l : lAbdo et al.ll2010d : iD'Ammando et al.ll201ll ). In addition, the appar- 
ent correlations between some major 7-ray flares with the coherent rotations of the optical 
polarization vector (swings by > 360°), as well as with the structural changes i n the inner ra- 
dio jet, have been discussed by lMarscher et al.l fl2010l ) and lOrienti et al.l (120 121 ) as important 
findings supporting the 'far-dissipation zone' scenario for the source. It is not clear, how- 
ever, if this model can be easily reconciled with the short variability timescales found in our 
analysis for the powerful 7-ray outbursts dissipating the bulk of the total kinetic luminosity 
of the jet in the s ystem. A substantial sub-structure of the outflow, consisting of a 'jets-in- 
jet' configuration (iGiannios et al.ll2009l ). or of a highly turbulent relativistic plasma passing 
through a standing shock (in which c ase one could e xpect an enhanced energy dissipation 
only in a small fraction of the jet; see Marscher 2012), may be possible solutionsEI 
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Table 1: Major 7-ray flares of PKS 1510-089. 



MJD 


^>100 McV 




(1) 


(2) 


(3) 


55853.5-55854.5 


14.86 ±0.89 


1.97 ±0.04 


55867-55869 


10.95 ±0.57 


2.21 ±0.04 


55872-55874 


8.39 ±0.44 


2.19 ±0.04 



(1) Dates of the three major 7-ray flux maxima m the daily-binned light curve of PKS 1510—089; (2) photon 
fluxes measured at the flux maxima in the units of [10^^ plicm^^ s^"'^], averaged over the specified time 
intervals; (3) the corresponding photon indices. 
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Fig. 1.— Daily 7-ray light curve of PKS 1510-089 during the period MJD 55834-55903 
analyzed in this paper. 95% flux upper limits are represented by triangles. Horizontal lines 
separating the three major flares are chosen arbitrarily just to guide the eye. 
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Fig. 2. — Fermi-LAT light curves of PKS 1510—089 around the time of the first major 7- 
ray outburst, binned in the intervals of 12 h, 6h, and 3h (upper, middle, and lower panels, 
respectively). 
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Fig. 3.- 



Same as in Figure[2]but for the second major 7-ray outburst. 
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Fig. 4.— 



Same as in Figure[2]but for the third major 7-ray outburst. 



